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Abstract 



d We analyze the spatial and energy dependence of the local density of states in a 



SNS junction. We model our system as a one-dimensional tight-binding chain which 



^ we solve exactly by numerical diagonalization. We calculate the dependence of the 

rH Andreev bound states on position, phase difference, gate voltage, and coupling with the 

• superconducting leads. Our results confirm the physics predicted by certain analytical 

C^ approximations, but reveal a much richer set of phenomena beyond the grasp of these 

H approximations, such as the metamorphosis of the discrete states of the normal link (the 

,JL normal bound states) into Andreev bound states as the leads become superconducting. 

o 

lHj 1 Introduction 

►^ The Josephson, or superconductor-norinal-superconductor (SNS) junctions, reveal some of 

00 the most interesting aspects of superconductivity. They are made of two superconducting 

^ (SC) layers and a normal link, and their physics is dominated by the SC proximity effect by 

fSj which a gap is induced in the spectrum of the normal link. When taken out of equilibrium, 

1^ a supercurrent (Josephson current) flows through such a junction p. The appearance of this 

O current seems paradoxical, as no current is expected to flow trough a gapped normal link. To 

^^ solve this, it has been proposed that discrete states, called Andreev bound states (ABS), are 

formed inside the gap [2]. The supercurrent is carried by these states, much like a regular 
current is carried by the Landauer-Buttiker transmission channels in a ballistic system; thus, 
the Josephson current is proportional to the number of ABS. This current is also proportional 
to the phase derivative of the energy of the ABS states: the larger the sensitivity of an ABS 
to the phase difference between the superconductors, the larger the Josephson current carried 
by it. While these states have been predicted quite a long time ago [2], they had only been 
measured indirectly via transport experiments [31 S], until their recent direct detection via 
spectroscopic techniques [SI E] . 

The ABS have been explored theoretically via the Bogoliubov-de-Gennes (BdG) equations, 
the Andreev approximation [71|8], as well as other approaches [9j. These techniques give the 
ABS spectrum in various limits, but cannot capture all the complex effects associated with 
the flnite length of the link, the eventually asymmetric, imperfect coupling between the link 



and the superconductors, the gate voltage, etc. Moreover, the dependence of the density of 
states in the hnk on position has been addressed only in a few works, mostly in relation to 
the form of the ABS wave-functions [TU] . 

In this work we calculate the full ABS spectrum, including the limits in which the analytical 
approximations described above do not work, by using a lattice one-dimensional tight-binding 
model [11] that we diagonalize numerically. This procedure allows us to obtain the dependence 
of the ABS spectrum on position, as well as on other parameters such as the SC order param- 
eter or the coupling with the leads. By tuning these parameters continuously we capture the 
transition between the SC state and the normal state and the way this is manifested in the 
ABS spectrum. 

Furthermore, we elucidate the relation between the normal state bound states of the link 
(NBS) and the ABS. The NBS spectrum is made of discrete sharp levels, whose spacing is 
controlled by the length of the link; when the normal link is coupled to metallic leads the NBS 
levels become wider, and the spectrum becomes energy independent for a perfectly-coupled 
system. Our most important observation is that the ABS are connected to, and evolve from 
the NBS. Thus we find that the number of ABS is proportional to the number of NBS falling 
inside the proximity-induced gap: the larger the ratio between the SC gap and the NBS level 
spacing, the larger the number of ABS. Moreover, the spatial structure of the ABS and the 
NBS is the same, and there is a one-to-one correspondence between each ABS and NBS in 
the spectrum! We confirm this correspondence between the ABS and the NBS spectra by 
comparing our numerical results with the results of analytical approximations [71 E]. 

The energies of the ABS are in general different from the energies of the NBS, and the 
energy difference depends on the coupling to the leads and the SC gap. We find that the larger 
the difference in energy between the ABS and the NBS, the larger the phase variation of the 
ABS. When the SC gap is much larger than the coupling with the leads, the NBS and the 
ABS have the same energy, and the dependence of the ABS position on phase is very small. 

In section 2 we present the tight-binding Hamiltonian that we use to describe the system. 
In section 3 we present the dependence of the density of states (DOS) on position and energy 
for non-SC contacts, focusing in particular on the dependence on the coupling with the leads. 
In section 4 we analyze the SC state, focusing on the dependence of the DOS on the coupling 
with the leads and the SC gap (section 4.1), on the gate voltage (section 4.2), on the phase 
difference between the two SC (section 4.3), on position (section 4.4), and in section 4.5 
we describe the DOS in a double-dot configuration. In section 5 we discuss the Andreev 
approximation (section 5.1), the "large-gap" approximation (section 5.2) and the wide-band 
approximation (section 5.3), and compare their outcome to the numerical results obtained via 
the tight-binding model. 

2 Model 

We model our system as a one-dimensional chain of atoms with inter-atomic distance a, and 
with one electron per atom. We assume that the number of sites in each lead is S, while the 
number of sites in the link is A^. The electrons are allowed to hop between atoms, with a 
hopping parameters tij (see Fig. [I]). The number of electrons per atom can be modified by en 
electrostatic potential (gate voltage) which can depend on position, Vq. The superconductivity 



in the leads is modeled by an on-site pairing potential Aj [8]: 

H = ^ (Yl ^G^Lcia + Aid^di) - Y^ tijclcj + h.c, (1) 

i=l,2S+N o=t,4- <i,j> 

where 'Y2<i i> denotes summing over nearest neighbors. 

For simplicity we will assume that for i < S, Ai = A exp"**^/^, for S < i < S + N , Ai = 0, 
and for i > S + N, A^ = Aexp**^/^. Thus the superconducting order parameter is uniform 
in the leads, with a phase difference of (p between the two leads, and is zero inside the link. 
An entirely normal system is retrieved when we set A = 0. For most of the calculations we 
assume that in the leads the applied gate voltage is zero, while in the link the gate voltage is 
uniform, Vq = fi. 

We should note that we do not impose periodic boundary conditions on our chain, so that 
no confusion arises when incorporating into the model the superconducting gap and the phase 
difference between the two SCs. 

This Hamiltonian can be solved by writing it in a 2(25* + A^) x 2(25* + A^) matrix form (two 
states for each sites corresponding to the spin-up/-down) that can be diagonalized exactly 
numerically using techniques similar to those of [TT]. This diagonalization allows one to 
obtain the entire eigenvalue spectrum, as well as the corresponding eigenfunctions. Also, by 
inverting the {u + iS)l2{2S+N) — H matrix one can obtain numerically the retarded Green's 
function matrix, which in turns allows one to obtain the DOS as a function of energy and 
position. 

To compare the validity of this numerical method we compare the numerical results with 
the corresponding results for an infinite chain. Thus, for an infinite uniform system with 
hopping parameter t the Hamiltonian can be written as 

H= r \e{k) - f^]clc, (2) 

where e{k) = —2tcos{ka). The dependence of the DOS on energy for this simple model in the 
normal state can be calculated: 

piE) = -Im r " GRik, tu) = -Im T " ^^ oc ^^. ^^''"^'^ (3) 

As we can see in Fig.lllc), this form of the DOS is reproduced by the numerical calculations, 
even if the size of the considered chain is not too big (100 sites). We note that the bandwidth 
of an uniform infinite system with hopping parameter t is 4t, with allowed states of energy 
between (— 2t,2t). The form of the DOS is well reproduced numerically also in the SC state 
(seFig.gd). 

3 Normal state 

3.1 Single impurity 

We start our analysis by reviewing the simple example of an infinite chain with a single 
localized impurity. We model the impurity as an on-site delta- function potential, Vq. Friedel 
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Figure 1: a) One-dimensional chain, b) The energy dispersion for the infinite uniform system and 
(c) the corresponding density of states d)The density of states when the gap is non zero (A = 1). 
In c) and d) we have set t = 1, and we have performed the calculations for a uniform system of 100 
sites, with a roughly large (5 = 0.1 



oscillations are generated in the DOS by the scattering of the electrons by the impurity, as 
described in previous work [12]. Other types of impurities may be considered, such as a local 
modification of the hopping parameter, but the qualitative features presented here do not 
depend on the exact form of the impurity for energies and positions that are not too small. 
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Figure 2: The density of states as a function of position and energy in a one-dimensional chain with 
a single impurity. We have taken t = 1, Vq = 1, and 5 = 0.04. 



It is most instructive to plot the LDOS as a function of energy and position, such as in 
Fig. |2] The chosen energy range is (— 2t,0), below — 2t, the bottom of the energy band, the 
DOS goes to zero. This type of plot allows one to easily identify the Friedel oscillations, 
and study their dependence on energy. Thus, close to the bottom of the band, the period of 
the Friedel oscillations is very large, since they arise from scattering of electrons that have 
similar momenta, as can be seen from Fig. [T| At energies close to zero, the allowed momentum 
transfer in a scattering process is close to 2ko, where ko = 7r/2a, yielding Friedel oscillations 
with a period of 2a (double the atomic lattice constant), and alternating contrast in the DOS 
on neighboring atoms. 

3.2 Double impurity - quantum dot connected to leads 

We focus next on a one-dimensional chain with two impurities located at the sites S, and S+N. 
This can be viewed as a simplified model of a finite-size system connected to metallic leads 
via less-than-perfect contacts. One can describe the impurities either by on-site potentials, as 



we discuss in section 3.3 or by local modifications of the hopping parameter, both yielding 
similar qualitative results. In what follows we model the impurities as modifications of the 
hopping parameter. The bad-coupling regime is achieved in this setup when the hopping 
between the normal link and the leads goes to zero. On the other hand, the good coupling 
regime is not achieved when the hopping between the link and the leads tj^p is infinite, but 
when it becomes equal to y/tieadstdot, where tieads and t^ot are the bulk hopping parameters in 
the leads, and in the dot respectively pj. If tdot = Uead = t (the dot and the leads are made of 



the same material), this prescription yields a perfect coupling for timp = t (uniform system), 
as exDected. 
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Figure 3: The DOS as a function of position and energy in a one-dimensional chain with two 
impurities. We have taken S = 50, N = 20, 6 = 0.05, t' = 0.2, and t = 1. 



We can see in Fig. [3] that the DOS exhibit peaks which correspond to forming of bound 
states in the finite-size system, and to the quantization of the allowed momenta k = rni/L. 
Consequently, the allowed energy values are En = —2tcos{n7ia/L) = —2tcos[7m/{N + 1)] , 
where N is the number of sites. Ignoring the spin, there are N bound states {n = 1, ...A^) for 
the entire bandwidth (— 2t, 2t). Note that if the number of sites is odd, the spectrum contains 
a zero-energy mode, (obtained when n = N + 1/2), however, if the number of sites is even, 
the minimum of the absolute value of the energy is |£'Ar+i/2| = 2sin{7r/[2(A^ + 1)]}. This 
dependence of the DOS on position and energy has been studied before in one-dimensional 
systems, both experimentally and theoretically, for both interacting, and even non-interacting 
systems [13]. 

The particular model that we are using allows us to study easily the transition between the 
Fabry-Perot regime (the central region is coupled almost perfectly to the leads) and the bad- 
coupling "Coulomb blockade" regime [H]. However, we should note that we cannot capture 
here the actual Coulomb blockade phenomenon, nor the extra separation between energy levels 
induced by the Coulomb interactions, since such interactions are not included in our model . 
To study the transition between the good-coupling and the bad-coupling regimes we plot the 
density of state for a site on the link as a function of energy for various values of the coupling 
with the leads. In Fig. |4] we note that the sharp peaks observed in the bad-coupling regime 
evolve into wider ones with increasing the coupling, such that the DOS has a quasi-sinusoidal 
energy dependence in the well-coupled regime. 
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Figure 4: a) Two-dimensional plot of the DOS on the site S + N ji as a function of energy and 
coupling with the leads b) The DOS on one central site as a function of energy for three different 
values of the coupling with the leads: t' = (full line), i! = 0.5 (dashed line), and t' = 0.8 (large- 
dashed line). The hopping parameter in the bulk was taken to be t = 1. The other parameters are 
b = 0.05, S = 100, and N = 20. 



3.3 Comparison with the continuous model for an infinite system 

The single- and double-impurity problems can be solved also using the T-matrix approximation 
when the length of the wire is infinite. For simplicity we consider a one-dimensional system 
with a quadratic energy dispersion: 

H= n e(A;)4c,, (4) 

J —Tr/a 

where e{k) oc /c^. The quadratic dispersion is a good approximation to the energy dispersion 
of the tight-binding chain for Fermi energies close to the bottom of the band (almost zero 
filling). Note that in order to model the tight-binding chain at half filling we need to study a 
system with linear dispersion whose spectrum includes right-moving and left-moving particles 



/n/a 
e(^)(ckcRfc-c[fcCLfe) (5) 

-n/a 



where c{x) = e^^''^CR{x) + e~*'^^^CL(a;) and cru and cik are the respective Fourier transforms 
of the fermionic operators. The retarded Green's functions can be written for this system as 
a function of energy and position: 

Gr{x, uj)= f dk^-—e'^^ (6) 

J CO + id — k^ 

when we work in units in which the mass of the quasiparticle is equal to 1. Using this form 
for the Green's functions, and assuming that the impurity potential is assumed to be a delta- 
function on-site potential V{x) = V6{x)c^{x)c{x), the LDOS dependence on position and 
energy can be determined straightforwardly using the T-matrix approximation |15l : 



p{x,oj) = -lmGR{x,uj)GR{-x,uj)T{uj) (7) 

with 

^(-) - 1-J«(0..) («' 

In Fig. [5] we plot the resulting LDOS as function of energy and position. Indeed, As expected 
for a system with quadratric dispersion, at low energy we observe long-wavelength Friedel 
oscillations, while, at high energy, the DOS exhibits features qualitatively similar to those 
arising in a tight-binding chain at half filling, such as high frequency, short-wavelength {kp'^/2) 
Friedel oscillations. We can see that the numerical and analytical result retrieve qualitatively 
the same form for the Friedel oscillations for energies and positions that are not too small (at 
low energy the analytical calculation has problems converging, yielding the white cut-out of 
the LDOS figure). 

While not shown here, a similarly good agreement was obtained for the DOS in the presence 
of two impurities (V"(x) = V6{x)c^{x)c{x) + VS{x — L)c^{x)c{x)) using the T-matrix formalism 
in Ref. |16|. 
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Figure 5: The DOS as a function of position and energy in a one-dimensional chain with a single 
impurity. The impurity potential is set to 1/ = 1. 



3.4 Double quantum dots 

We can easily adapt our formalism to study double-quantum-dot systems by modifying the 
hopping between the sites S + N/2 and S' + A^/2 + l. In Fig. [6] we are plotting the DOS on one 
of the dots as a function of the two gate voltages applied separately on the two dots. If the 
two dots are completely decoupled, this plot yields a grid of horizontal lines - the states of one 
dot do not depend on the gate voltage applied on the second dot. However, when the dots are 
coupled, "transitions" occur between the horizontal lines, giving rise to rounded steps. These 
phenomena have been discussed in detail in the literature [T^, and we confirm here that a 
tight-binding model recovers a similar form for the dependence of the DOS on the two gate 
voltages. 



4 Superconducting state 

Our main interest is to investigate the formation of ABS in the normal link when the leads 
become superconducting. A vast literature exists detailing the properties of ABS, their as- 
sociated wavef unctions, the corresponding Josephson current, etc.|2l |5|, [6l [TJ [HI |9] . While the 
ABS have been detected experimentally indirectly in transport measurements [31 S], a direct 
spectroscopic measurement has not been done until recently [HI [6]. Most of the ABS theo- 
retical studies have been performed using the Bogoliubov-de-Gennes equations and different 
approximations [3 El E] . Here we use a tight-binding approach, which is more versatile and 
allows us to modify easily any parameter in the system (coupling with the leads, be it sym- 
metric or asymmetric, phase difference, gate voltage, position), thus giving access to a much 
richer picture than the one obtained analytically. Moreover, our technique permits to see how 
the ABS evolve from the bound states existing already in the non-superconducting system, 
and to study the spatial dependence of the ABS. Last, but not least, it provides an easy tool 
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Figure 6: The DOS on one dot as a function of the two gate voltages on the two dots. The parameters 
are t = 1, Unter-dot = 0.3, t' = 0.5. 



to take into account the effects of disorder. 



4.1 The dependence of the ABS on the coupUng with the leads, as 
well as on the superconducting gap 

We begin by considering a one-dimensional chain with two impurities (similar to the one 
discussed in the previous section) for which the superconducting order parameter in the leads 
is non-zero. As described in previous work, a superconducting gap is induced in the spectrum 
of the normal link, and discrete ABS are formed inside this gap. In Fig. [7]we plot the spectrum 
of a system with A^ = 20 and S = 200 for an isolated system (zero coupling with the leads) 
(dotted line), as well as for a coupling of t' = 0.5 (full line). When the coupling is non-zero, 
the spectrum of the system is modified: gap-edge-like features appear at energies close to 
A = 0.5, and the position and shape of the NBS peaks change. Thus the peaks inside the gap 
change position but remain sharp - these peaks correspond to the ABS. On the other hand, 
the peaks outside the gap are shifted by much smaller amounts, but are smoothed out by the 
coupling with the leads, same as in the fully normal state. While the energies of the ABS 
are different from those of the NBS, there seems to be a one-to-one correspondence between 
each NBS and the ABS. To understand this issue better, in Fig. |8] we plot the DOS as a 
function of energy on the x-axis, and, on the y-axis, the coupling with the leads (top panel), 
and the superconducting gap (bottom panel). Fig. [s] allows us to understand the evolution of 
the spectrum between the normal state and the superconducting state. Thus in the bottom 
panel we note that in the normal state (A = 0) there are roughly four allowed bound states 
in the analyzed energy interval, which are broadened by the coupling with the leads. When 
the value of the gap becomes comparable with the energy of the first state, this state "enters" 
the gap, becoming an ABS. When the gap is increased even more, the second state also enters 
the gap, etc. For small gaps the ABS are pretty much "stuck" to the gap edge, while when 
the gap increases their energy evolves asymptotically towards a given value. By comparing 
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Figure 7: The DOS as a function of energy for a system witli S = 200, N = 20, and 5 = 0.02. The 
site for which the DOS is evaluated is S" + N/2 = 210. The hopping term t has been set to 1. The 
coupling with the leads is set to zero for the dotted line, and to t' = 0.5 for the full line. The SC 
gap is set to A = 0.5. 



the lower panel with the upper panel, which describes the DOS as a function of energy and 
the coupling with the leads, we can see that the positions of the ABS at asymptotically large 
gap are the same as the positions of the NBS (the NBS spectrum corresponds to t' = in the 
upper panel, and the NBS positions are denoted by the black dashed lines). 

Fig. |8] - upper panel - describes how the ABS spectrum is modified when the value of 
the coupling with the leads is modified. The states outside the gap (the position of the gap 
A = 0.8 is marked in the upper panel by the white dotted line) become wider when the 
gap increases, while the states inside the gap have their position slightly modified by the SC 
coupling, but their widths are unchanged. 

From our results one can thus establish a direct connection between the NBS and the ABS 
spectra. For example, if there is no NBS in the energy window delimited by the SC gap, no 
ABS will form, or the ABS will be stuck to the gap edge. Moreover, the number of ABS is 
directly related to the number of NBS in this energy window. We will discuss in more detail 
in section[5]the analytical interpretation of the correspondence between the NBS and the ABS 
and in what limits this correspondence can be established. 



4.2 ABS dependence on gate voltage 

One of the most interesting aspects of the ABS is their dependence on gate voltage. It is this 
dependence that allowed one to unambiguously identify the ABS experimentally O [6]. In 
Fig. [9] a) we plot the DOS as a function of energy and gate voltage for a system of S = 100 
and A^ = 20, with t = 1, t' = 0.5, and a SC gap of A = 0.2. For this system the spacing 
between the normal bound states is El = At/N = 0.2. We note that at yU = (along the 
dashed line in the figure), two ABS form, one at positive energy and the other at negative 
energy, consistent with the fact that the gap is of the order of El. Away from /i = the 
states split each into two ABS, which evolve forming bell-shaped/sinusoidal features, and the 
positive-negative energy symmetry is broken. Outside the gap the position of the states is 
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Figure 8: The DOS as a function of energy on the x-axis, and, on the y-axis, the couphng to the 
leads {t') (top panel) and the superconducting gap (A) (bottom panel). The coupling with the leads 
is assumed to be symmetric. We have taken A = 0.8 for the upper plot, and t' = 0.5 for the lower 
plot, S = 150, N = 20, and 6 = 0.02. The site for which the DOS is evaluated is S + N/2 = 160. 
The hopping term t has been set to 1. The positions of the NBS are given by the black dashed lines, 
while the gap is indicated in the upper panel by the white dotted line. 
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linear with gate voltage, same as in the normal state; these normal states transform into ABS 
when they touch the gap edge. Similar features have been described using a simpler model in 
0. 
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Figure 9: The DOS as a function of energy and gate voltage. We have set a) A^ = 20, S* = 100, 
t = 1, t' = 0.5, A = 0.2, and 6 = 0.02; b) same as in a) plus a phase difference between the two SC of 
27r/3; c) same as a) plus an asymmetric coupling of the link with the two SCs, t^ = 0.25, t^ = 0.75; 
d) same parameters as in a) except for A = 0.5. The white dashes lines A and B correspond to a 
gate voltage of /u = and respectively fi = n/N. 



We find that, consistent with previous observations, when the value of the gap increases, 
the number of ABS also increases. Also when the phase difference between the SCs is non- 
zero, or in the presence of an asymmetric coupling with the two SCs, the ABS at /i = split. 
The ABS also split when the gate voltage is non-zero. We should also note that we do not 
observe any zero-energy crossing of ABS levels. As described in [5], such crossings, which 
have been observed experimentally, are the effect of Coulomb interactions and of the lifting 
of degeneracy between the spin-up and spin-down energy levels. In our tight-binding model 
such effects are not included and indeed, the crossings do not appear. 
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4.3 ABS dependence on phase 

We now turn to the study of the dependence of the ABS on the phase difference between the 



SCs. In Fig. 10 a), b) and c) we plot the DOS as a function of energy and phase difference, 
for a A = 0.5 and /x = 0, corresponding to the phase dependence of the A spectrum in Fig. |9] 
We can see that the two degenerate = ABS spht for non-zero (f). The positions of the ABS 
vary sinusoidally with 0, the largest amplitude of this variation is achieved when the coupling 
with the leads is perfect t' = 1 (a). Moreover, consistent with observations in the previous 
section, the degeneracy of the two ABS at = is lifted for an asymmetric coupling with the 
leads t'^ = 0.25, t'j^ = 0.75. We should also note that even some of the states outside the gap 
and close to the gap edge have a non-zero phase dependence, weaker than that of the states 
inside the gap. 

For the generic parameters we have considered, the dependence of the ABS on flux changes 
qualitatively when the parity of the number of sites in the normal link changes. The difference 
stems from the effects of this parity on the NBS spectrum. As described in the previous section, 
if the number of sites is odd, a NBS zero energy state is formed, while if this number is even, 
there is no zero-energy NBS. This is equivalent to saying that the even and the odd spectra 
are shifted with respect to the zero-energy level by a gate voltage fi = n/^N + 1) equal to the 
difference in energy between the lowest energy eigenstates of the odd and even systems. For 
a large enough A^, we can recover an odd spectrum from an even spectrum by a simple gate 
voltage shift. The spectrum in Fig. [lo] (d) evaluated for N = 20 at fx = n/N (corresponding 
to the B line in Fig. |9]) should thus be the same as the /x = spectrum of a system with 
N = 21. 

4.4 ABS dependence on position 



In Fig. 11 we plot the DOS as a function of energy and position for a) the normal state, and b) 
the SC state. We observe the SC gap in the leads, the proximity-induced gap in the link and 
the corresponding ABS. Our analysis reveals that there is a one-to-one corresponence between 
the NBS and the ABS, such that for each ABS there is a NBS with the same spatial structure. 
Thus, each ABS evolves in the presence of the SC leads from an NBS. The main difference 
between an ABS and its NBS counterpart is the phase dependence of the ABS which is more 
important for the ABS whose energies differs the most from the NBS energies. 



In Fig. 12 we zoom in on the lower-energy states: in a) the parameters are identical to the 



ones in Fig. 11 , in b) we have introduced a non-zero flux = 27r/3, in c) the coupling with the 
leads is asymmetric, t'^ = 0.25, t'j^ = 0.75, while in d) the gate voltage is non-zero fi = —0.05, 
A = 0.2 and 6 = .01. We note that the two states obtained from the splitting of an ABS 
when 7^ have the same spatial structure. However the two states obtained when the ABS 
split due to an asymmetric coupling have a different spatial structure. Furthermore, all the 
four states (negative and positive energies) arising at /x 7^ have the same spatial structure. 



as it can be seen in Fig. Hi). Last, but not least, our analysis confirms that the ABS "leak" 



into the normal leads, and that this leak is larger when the ABS are closer to the gap edge. 

4.5 ABS in double quantum dots 

By modifying the coupling between the sites S + A^/2 and S + A^/2 + 1 we can split the normal 



link into two quantum dots. In Fig. [13] we plot the dependence of the DOS on one dot as a 
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Figure 10: The DOS as a function of energy and phase difference on the site S+N/2. The parameters 
are S = 100, N = 20, t = I, 6 = 0.02, A = 0.5, and a) t' = 1, b) t' = 0.5, c) t'^ = 0.25, t'^ = 0.75 
d)sanie parameters as a) but for non-zero gate voltage /i ~ tt/{N + 1). The spectra in a), b) and c) 



correspond to the hne A in Fig. M and the spectrum in d) corresponds to the hue B in Fig. 10 
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Figure 11: The DOS as a function of energy and position for a) t = 1, t' 
N = 20, S = 100, = 0, /x = 0; b)sanie parameters except for A = 0.5. 



0.5, A = 0, 5 = 0.04, 
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Figure 12: The DOS as a function of energy and position. In a) the parameters are identical to the 



ones in Fig. 11 , in b) we have introduced a non-zero flux (p = 27r/3, in c) the coupling with the leads 



is asymmetric, t'^ = 0.25, i^ = 0.75, while in d) the gate voltage is non-zero /i = —0.05, A = 0.2 and 
5 = .01. 
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function of the two gate voltages on the two dots. We can see that the picture is very similar 
to the normal state picture in Fig. |6} 




-1.0 

-1.0 -0.5 



0.0 0.5 1.0 



Gate voltage 1 

Figure 13: The DOS on one dot as a function of the two gate voltages on the two dots. The 
parameters are set to t = 1, t' = 0.5, titter-dot = 0.5, A = 0.5, and 5 = 0.04. 



5 Discussion 

To understand better our results we discuss three particular limits in which the BdG equations 
can be solved analytically to obtain the energies of the ABS, and we compare these analytical 
results to our tight-binding results. 



5.1 Low-energy Andreev approximation 

If the lowest-energy ABS form at energies much smaller than the gap, t/N ^ A < t, one is in 
the regime of validity of the Andreev approximation. The ABS energies (for perfect Andreev 
reflection) can be obtained in this limit from the quantization of the phase acquired inside the 
link [7]: 

2 arccos(E/A) - {q^ - qh)L ± = 2mx (9) 

where Qe and qu are the corresponding momenta for the electron and the hole with energies E 
and respectively —E, where E = —2t cos[{kp + qe)a] = 2tcos[(/ci? — qh)o]- 

For energies much smaller than the gap {E ^ A, i.e. arccos(-E'/A) ^ 7r/2), this becomes: 

(ge-gh)^ = vr(2n-l)±0. (10) 

Given the band structure presented in Fig. 1, at energies much smaller than t, and for chemical 
potentials ^ <t^ t, the band can be linearized around the Fermi energy, with kp ~ 7r/2a + 
fj,/{2ta), and qe — Qh ^ —E /{at) independent of /i. We thus have 



E 



-t[7r(2n-l)±0]/(A^ + l) 



(11) 
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where we have used the fact that L = {N + l)a. We can compare this quantization relation 
to the one used to obtain the NBS: 

2keL = 2Tin (12) 

where E = fi — 2t cos{2kf.a) , which yields, as discussed in Section 2.1, 

E = fi-2tcos(^j^^ (13) 

The energies of the NBS are different if N is even or odd. Thus, at low energy [E ^ t) 

Eodd = /i - 2n7r/{N + 1) (14) 

and 

^e.en = /i-(2n + l)vr/(iV + l). (15) 

Hence, for A^ even the NBS spectrum is identical to the ABS spectrum evaluated at = 0, 
while for A^ odd, the spectrum of the NBS is the same as the spectrum of ABS at = tt. This 



is confirmed numerically; in Fig. 14 we plot the ABS dependence on for N = 100, t = t' = 1, 
and A = 1. Indicatively, in dashed lines we give the NBS energies for A^ = 100 and in dotted 
lines the NBS energies for N = 101. 




0.00 

0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Figure 14: Dependence of the DOS on the site 5" + 1 on energy and the phase difference between 
the SCs. The parameters A^ = 100, S = 20, A = 0.3, t = t' = 1, and 5 = 0.01 are chosen so that 
we are in the Andreev Umit. The dashed Hnes and the dotted Unes describe the NBS spectrum for 
A^ = 100 and A^ = 101 respectively. 



5.2 Large-gap approximation 

In what follows we analyze the limit in which the gap is much larger than the hopping pa- 
rameter t. This limit is quite unphysical, but it is interesting to study it analytically, since it 
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reveals the mechanism via which the ABS can be identified with the NBS for a specific range 
of parameters. We begin this analysis by a numerical study, focusing on the lowest-energy 
ABS, whose energy is much smaller than t. When t and A are comparable, the system is 
described quite well by the formalism presented in the previous section. A typical dependence 



on phase in this limit has been presented in Fig. [14} When A is increased above t, we can 

100, the amplitude of the variation of the ABS with the phase is 
= TT. When the gap becomes even larger, the amplitude of 



see in Fig. [T5| that, for N 
decreased, and a gap opens at 
the oscillations goes to zero, the ABS become phase-independent and their energies converge 
towards the NBS energies. While not depicted here, for A^ odd, the ABS energies also con- 
verge at large gap towards the energies of the NBS, which, when A^ is odd, are identical to 
the small-gap ABS energies at = vr. 




0.00 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 



Figure 15: Dependence on the DOS on energy and phase for a A^ = 100 link, when the system is no 



longer in the Andreev limit. This is achieved by using the same parameters as the ones in Fig. 14 
and increasing the gap to A = 2. We can see that the amplitude of the variation with the phase is 
reduced, and gaps open at (/> = vr. Same as above, the dashed lines and the dotted lines describe the 
NBS spectrum for A^ = 100 and N = 101 respectively. 



We can understand the origin of this phenomenon by finding the ABS energies exactly 
for an extra-small system with one SC site in each lead, and one normal site. When an 
expansion of the energies of the ABS is performed at large gap to first order in t/A, we find 
the allowed energies to be independent of the phase difference, and their values to be given 
by £^ = ±/x, ±A. Indeed, the E = n is the energy of the NBS for a single site at chemical 
potential fi. The extra energy in the spectrum, E = —fj,, has a negligible weight, so it does 
not give rise to a peak in the DOS. 



5.3 Wide-band approximation 

We can also work in the limit of a; ^ A ^ t. In this limit one can use the wide-band approxi- 
mation [T8| IT9] which allows one to take into account the superconducting leads by a adding to 
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the Green's function matrix of the normal hnk a self-energy resulting from integrating out the 
leads. In the limit u <^ A, the formalism presented in [181 [IS] yields self energy contributions 
for the right and left leads in the form of 2 x 2 matrices in the BdG basis, whose only non-zero 
elements are off-diagonal and are given by iF/j/^.e"'^^''^ . The DOS can be obtained exactly in 
this limit for a link made of one site or two sites. For a one-site link we find: 

E = ±T cos(0) (16) 

while for a two-site problem 

E = ±v/F2 + t2±2Ftsin(0/2) (17) 

We can see that in this limit the energy of the ABS is dominated by the coupling with the 
leads, rather than by the intrinsic energies of the NBS. 

6 Conclusion 

We have analyzed the spectrum of an SNS junction, focusing in particular on the dependence 
of the DOS on energy and position. We have used a tight-binding numerical method that has 
allowed us to study a broad range of parameters, and thus to go beyond previous analytical 
approaches. Our calculations yield similar features to those recently measured in SNS junc- 
tions with carbon nanotubes |5j and graphene quantum dots [6], which further confirms that 
the generic parameters used in our numerical calculations are realistic. Since these parameter 
values do not correspond to any limit accessible via analytical approximations, we believe that 
the full physics of the ABS in such a system can only be obtained numerically. 

Furthermore, our analysis reveals a one-to-one correspondence between the NBS states 
and the ABS states, whose spectra become identical when the the SC gap is much larger than 
the bandwidth of the normal link. In such a limit one may hope to obtain access to the NBS 
spectrum (which is in practice masked by the coupling with the leads) via a measure of the 
ABS spectrum. 

While not presented here, we have also verified that our numerical tight-binding results 
are consistent with results obtained by solving the corresponding BdG equations numerically. 
The analysis required to obtain such solutions is however much heavier, and less elegant, and 
we believe that the tight-binding method presented here is faster and more versatile. It would 
be interesting to generalize our calculation to include the effects of disorder, and to test the 
stability of the ABS for a disordered system. 

Last, but not least, our analysis provides a criterion for the spectroscopic observation of 
the ABS in a given clean ballistic system, based on the energy spacing of the NBS. Thus, if 
the distance between the NBS levels is much smaller than the SC gap, the ABS may be too 
close to be resolved in energy. Reversely, if the distance between the NBS is very large, the 
ABS will tend to stick to the gap edge, except when the NBS approach the Fermi level and 
one pair of ABS at a time will be resolved inside the SC gap. 
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